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Thermal gravity, black holes and cosmological entropy 

Stephen D.H. Hsi:Q and Brian M. Murraj0 
Institute of Theoretical Science, University of Oregon, Eugene OR 94703-5203 

Taking seriously the interpretation of black hole entropy as the logarithm of the number of mi¬ 
crostates, we argue that thermal gravitons may undergo a phase transition to a kind of black hole 
condensate. The phase transition proceeds via nucleation of black holes at a rate governed by a 
saddlepoint configuration whose free energy is of order the inverse temperature in Planck units. 
Whether the universe remains in a low entropy state as opposed to the high entropy black hole 
condensate depends sensitively on its thermal history. Our results may clarify an old observation of 
Penrose regarding the very low entropy state of the universe. 


I. INTRODUCTION 

Years ago Penrose noticed that the universe must have 
begun in a very low entropy state Q. By considering 
the entropy of black holes, he argued that the current 
state of the universe has significantly lower entropy than 
the maximum possible entropy state. For example, while 
holding the number of baryons fixed one could increase 
the total entropy tremendously by letting matter collapse 
into black holes Indeed, it seems that while the mat¬ 
ter degrees of freedom were born hot, i.e., in a maxi¬ 
mum entropy thermal state, the gravitational degrees of 
freedom were born in a very special low entropy state. 
Interpreting entropy as the logarithm of phase space vol¬ 
ume, a low entropy state is an exponentially unlikely 
state and hence can only result from fine-tuned initial 
conditions 0. Reasoning along these lines suggests that 
spacetimes with numerous horizons, perhaps resembling 
a dense agglomeration of black holes, occupy an expo¬ 
nentially larger fraction of gravitational phase space than 
smooth spacetimes like the usual Friedmann-Robertson- 
Walker TRW) cosmologies. For related discussions, see, 
e-g-, |l,p and references contained therein. 

One may ask whether special initial conditions at the 
Planck scale are sufficient to produce the low entropy 
universe we see today. It might be the case that interac¬ 
tions with thermal matter in the early universe inevitably 
cause the gravitational degrees of freedom to thermalize 
as well. Such a thermal state, assuming ergodicity of 
gravity, would likely evolve to a configuration of much 
higher entropy, and hence a cosmology very different from 
the one we observe. 

Since black holes are our only hint at the highly en- 
tropic configurations of gravity 0, they should play a 
prominent role in the transition from low entropy to high 
entropy spacetimes. In this paper we suggest a spe¬ 
cific mechanism involving the nucleation of black holes 
from a thermal graviton state. We note that the cor¬ 
responding nucleation rate from a thermal matter state 
is much smaller, and probably irrelevant cosmologically. 
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The mechanism we describe provides a plausible means 
by which Penrose’s ergodic evolution could proceed. We 
examine whether the transition to a new, highly entropic, 
phase of condensed black holes can occur in standard big 
bang cosmology. The result depends sensitively on the 
thermal history of the universe at early times. Moreover, 
the relevant energy scales are all higher than the energy 
scale at which an inflationary epoch is usually assumed to 
take place. Therefore, we are considering a phase tran¬ 
sition which only may take place before and not after 
inflation. Presumably, the probability for a given patch 
to inflate would be affected by whether or not that patch 
has undergone a phase transition to the high entropy 
phase. 

We should note that there is not a consensus on the is¬ 
sue of whether gravity is ergodic nor on the interpretation 
of the gap between the maximum allowed and the actual 
entropy in an FRW spacetime. Tipler showed that 
under a reasonable set of assumptions, closed universes 
are technically not ergodic, i.e., there is no Poincare re¬ 
currence. Moreover, Barrow 0 has pointed out that in 
a spacetime restricted to be FRW there is necessarily 
an entropy gap, i.e., the entropy in thermal radiation is 
much less than the entropy associated with a black hole 
of horizon size. However, the phenomena we described 
in the previous paragraph, which are investigated in this 
paper, are independent of these larger questions about 
general relativity. That is, the mechanism by which black 
holes are nucleated occurs on sub-horizon time and length 
scales. The statistical approach we take below is justified 
by the presence of a thermal bath of gravitons or other 
particles, whose existence is not in dispute. In this sense 
we do not require any assumption of ergodicity, except in 
some small sub-horizon patch. 

In Sec. ITTIwe consider black hole nucleation in a system 
of thermal gravitons and compare to a thermal system of 
matter. In Sec. cni we show that gravitons may ther¬ 
malize in the early universe even if they started out cold. 
We determine the conditions necessary for a phase tran¬ 
sition to a black hole condensate via percolation in Sec. 
Cvl Finally, in Sec. m we relate these results to Pen¬ 
rose’s observation. We use Planck units throughout, i.e., 
h = c = G = kB = ^. 



2 


II. STATISTICAL MECHANICS OF 
GRAVITONS 

Consider a box of hot gravitons. The probability for a 
fluctuation to lead to a black hole of radius R is 

P{R) ~ (1) 

where E = R/2 is the energy of the black hole, and 
N is the multiplicity of microstates which, when coarse 
grained, appear as a black hole of radius R. The proba¬ 
bility can be written as 

P{R) ^ (2) 

where F = E — TS is the free energy, and S = In V. We 
assume that black hole entropy is accounted for by grav¬ 
itational microstates, as suggested by results from st ring 
theory M- Using the Bekenstein-Hawking formula |l2ll 
for black hole entropy ^bh = ^/4, where A = AttR'^ is 
the area, we see that 

F{R) = R/2 - tvTR^. (3) 

Strictly speaking, we want the free energy relative to 
that of hot, flat space. This means we should subtract 
from 0 a correction Fq{R) ~ —R?T'^. It is easy to see 
that, near the saddlepoint found below, Fq is a negligible 
correction as long as the saddlepoint radius i?* is much 
smaller than the horizon size. 

The radius that maximizes the free energy, the sad¬ 
dlepoint radius, is given by R„ = . We obtain 

1 i?* <C where E[~^ ^ l/T^ is the horizon size 

for a radiation dominated FRW universe, so our anal¬ 
ysis encounters no difficulties from quantum gravity or 
causality. The corresponding maximum free energy is 
(see Fig. 0 

F* = (167rT)-^ (4) 

At the saddlepoint radius the black hole temperature 
is just equal to that of the heat bath. Black holes with 
0 < F < i?* shrink to zero size, leaving a weak-field 
phase with a thermal population of graviton states (grav¬ 
ity waves) on a smooth background metric. However, 
for F > F*, black holes grow without bound (they are 
colder than the environment), and for F > Fq = {2'kT)~^ 
the free energy is actually negative, less than that of 
F = 0. This instability may indicate a new nonpertur- 
bative phase of gravity, which is not asymptotically flat, 
and in which spacetime is all or partially filled with black 
holes. Such a phase is highly entropic and occupies an 
exponentially larger phase space volume than the smooth 
weak-field phase. 

The nucleation rate for supercritical black holes (which 
might be thought of as bubbles of the new nonperturba- 
tive gravitational phase) is controlled by the free energy 
at the saddlepoint, as in the usual case of a first order 
phase transition (for early papers on nucleation theory, 
see ^3)- This yields 


F(R) 



FIG. 1: Free energy versus black hole radius for a gas of hot 
gravitons. Black holes of size F < F* shrink, leading to a 
weak-field gravity phase. For F > F*, however, black holes 
grow, possibly giving rise to a new nonperturbative phase of 
gravity. 


(Strictly speaking, the dimensional prefactor could be 
modified by subleading terms in the exponent.) The 
physics is similar to that of nonperturbative baryon 
number violation in the standard model. There, the 
rate is controlled by the free energy of the electroweak 
sphaleron, which is the saddlepoint confi gur ation sepa¬ 
rating vacua of different winding number |l4l . 

Now consider a box of hot photons. If a fluctuation of 
size F and energy E satisfies E > R (we ignore factors 
of order one), it will inevitably evolve into a black hole 
0. Therefore, Eqns. ( 0 ) and @ for the probability for a 
fluctuation to lead to a black hole of radius F still hold. 
In this case, however, the entropy is not proportional to 
A. In order to evaluate the multiplicity iV, we make use 
of a bound on the entropy of a region of size F filled with 
thermal radiation, originally derived by’t Hooft Eg. By 
noting that matter in thermal equilibrium, energy scales 
as F ^ F^T^, and entropy as S ^ R^T^, and further 
requiring that the system not have already undergone 
gravitational collapse, i.e., F < F, ’t Hooft obtained the 
bound S < (again, we ignore numerical factors). 

Matter configurations that lead to black holes saturate 
this bound. Therefore, we see that the free energy of 
relativistic matter configurations which evolve into black 
holes is of the form F(F) ^ F — TR?!"^. 

Once a fluctuation of sufficient size to lead to a black 
hole has occurred, the evolution is then governed by the 
same physics as in the original graviton case. We there¬ 
fore compare the multiplicities of configurations of pho¬ 
tons and gravitons, and iVg, respectively, that will 
lead to black holes of critical size F*. For temperatures 
well below the Planck scale, A ^ and so fluctua¬ 

tions of critical size in the photon gas are suppressed rel¬ 
ative to the case of a graviton gas by a factor of roughly 




exp 


I 


167rr2 


( 6 ) 


A(T) ~ 


(5) The key difference between the two cases considered is 
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the entropy limit on thermal degrees of freedom. Ordi¬ 
nary matter cannot achieve the entropic density of grav¬ 
itational degrees of freedom, under the assumption that 
black holes are coarse grained objects with gravi¬ 
tational microstates. Nucleation of black holes is much 
more likely if these gravitational microstates are ther¬ 
mally occupied (i.e., in a graviton heat bath), than if one 
starts with hot matter and cold gravitons. 


III. GRAVITON THERMALIZATION 


coupling between Xi and 7 . Again may be signifi¬ 
cantly below the Planck scale. 

Let us obtain a more careful estimate of in a 
specific model: long-wavelength quantum gravity (EFT) 
coupled to scalar QED (see, e.g., [13). For a tree level 
calculation, the relevant terms in the effective lagrangian 
are: 


C = \g\i - \g^''g>>^ 


( 7 ) 


Natural initial conditions for the universe might have 
both gravitons and matter in thermal equilibrium. In 
some cases, however, thermal matter leads to thermal 
graviton populations even when the gravitons are initially 
cold. We can investigate the thermalization of gravitons 
using a long- wavelength effective field theory (EFT) for 
quantum gravity llbl . If one is only interested in pro¬ 
cesses occurring at energies sufficiently below the Planck 
scale, as is the case here, there is no obstacle to using 
the standard EFT approach of including all terms in the 
effective lagrangian that are consistent with the symme¬ 
tries of the system, in this case general coordinate trans¬ 
formation invariance. 

Without knowledge of the fundamental theory of quan¬ 
tum gravity, we are not able to write down the Boltz¬ 
mann equation for the evolution of the phase space dis¬ 
tribution for the gravitational degrees of freedom. Boltz¬ 
mann heuristics, however, motivate using r(T) = H(T) 
as a reasonable criterion for freeze-out of a given particle 
species, in this case the graviton. r(T) is the interac¬ 
tion rate of gravitons with the heat bath at temperature 
T, and F[{T) is the Hubble expansion rate at tempera¬ 
ture T. Below the decoupling temperature Tdec the Hub¬ 
ble expansion rate is greater than the interaction rate of 
gravitons with the heat bath, and gravitons are decou¬ 
pled. 

Let us estimate the graviton decoupling temperature 
based on two different scattering processes. First, con¬ 
sider the interaction XiXi gg, where Xi is any one of 
N scalar, fermion or vector particles, and g is the gravi¬ 
ton. In a relativistic gas, Boltzmann heuristics suggest 
that the interaction rate is roughly TXiXi^gg ^ ntai, 
where rii and ai are the number density and cross sec¬ 
tion, respectively, of species i. At energies not too far 
below the Planck scale all particle species are relativistic 
so rii ~ . The matrix element for this process goes 

as and the Hubble rate as FI We see 

that the decoupling temperature obtained from consid¬ 
ering the process XiXi gg is Tdec For a 

model with a large number of particle species this scale 
may be well below the Planck scale. 

This estimate, however, does not take into account in¬ 
teractions between particle species. A similar argument 
based on considering the process XiXi 75 , where 7 
is another particle species, leads to a decoupling temper¬ 
ature of roughly ~ , where is the 


where = d^ + ieAfj_, F^i, = dfj^Ay — d^A^, and TZ is the 
scalar curvature. The metric is expanded about the 
Minkowski metric = diag(l, —1, —1, —1), i.e., 5 ^ 1 / = 
-b •\/327r/i^i,, and Feynman rules are obtained by the 
standard procedure. The cross section for the interaction 
(fxf) "fg is found to be 


327r 


( 8 ) 


where a = e^/47r. Assuming local thermal equilibrium, 
the scalar number density is = (2C(3)/7r^)T^. More¬ 
over, if there are Xg species of scalar particles, the gravi¬ 
ton decoupling temperature is 


^ dec 


5aW 


1/2- 


( 9 ) 


For consistency we must also estimate the temperature 
below which we can trust the above calculation, i.e., esti¬ 
mate the energy scale below which the EFT description 
is valid. Let us look at the next order contribution to the 
matrix element for this process: 


+ ■ ■ ■ ■ ( 10 ) 


Notice = C>(e-\/327rg), where q is some typ¬ 

ical energy scale of the interaction. In order to ob- 

C3'\ 

tain a conservative estimate of Alt,i , we will as- 

> 73 ’ 

sume that all diagrams interfere constructively. Then 
= 0{FLe{V3^q)^), where F is the number 
of graphs at this order and L is a loop factor associated 
with each graph. The EFT is valid for 

9 < 9pert = (327rFL)“i/^. (11) 

Taking F = 27 and L = l/ 47 r^ we get 9pert — 10“^. By 
comparing Eqs. 0 and (HU), we see that < qpert 
in the large Ng limit, due to the fact that ~ , 

while qpert is independent of Ng. Thus, there is a class 
of models (those with large numbers of matter fields) in 
which gravitons are copiously produced and interact fre¬ 
quently, at energy scales where the effective field theory 
applies. 

One could perform a similar calculation with a more 
realistic model for the matter content of the early uni¬ 
verse. At the energy scales of interest, say 10^® GeV, the 
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matter interactions may, for example, be described by a 
grand unified theory (GUT) m with Ng ^ 10^ and cou¬ 
plings not much smaller than unity. A calculation in such 
a model would likely yield the same conclusion, i.e., that 
there was an epoch in which even initially cold gravitons 
interacted strongly enough with the heat bath that they 
were themselves thermalized. 


IV. PERCOLATION OF BLACK HOLES 

In light of the previous sections, it is interesting to 
consider the possibility of a first order phase transition 
in an FRW universe with a thermal population of gravi¬ 
tons. Assuming that gravitons are in thermal equilibrium 
in the early universe between an initial time to and a 
later time ti, the volume fraction in the weak-field grav¬ 
ity phase is given by: 


p{to,ti) = exp 





In flat FRW spacetime 





( 12 ) 


(13) 


Here a is the scale factor, and b is the (constant) speed 
at which the black holes expand. Eq. is a standard 
formula from old inflationary cosmology 0 . b is gov¬ 
erned by relativistic physics and, therefore, should not 
be significantly smaller than unity. During a radiation 
dominated epoch the scale factor goes as a{t) ^ and 

temperature and time are related by t ~ O.Sg* _ 

The volume fraction (now as a function of temperature) 
is then of the form: 


p{To,Ti) = exp 



f{To,Ti) 


(14) 


We define the critical initial temperature Tq, so that if 
the initial graviton temperature is larger than Tg, the 
volume fraction in the weak-field phase is nearly zero. 
See Fig. [3 for a plot of p{To,Ti) for b = 10“^, 10“^, 
and 10“^, with Ti = 10“® and g* = 10^. Decreasing the 
expansion speed by two orders of magnitude has the effect 
of increasing T§ by less than a factor of two. Increasing 
the number of effective degrees of freedom or the cut-off 
temperature Ti both have similar effects on Tg, i.e., they 
cause it to increase by a relatively small amount. 

The critical initial temperature T§ is of order 10“^ to 
10“^, which from Eq. HU is the same order of magni¬ 
tude as the temperature at which quantum gravity be¬ 
comes perturbative. Thus, taking the initial temperature 
of the big bang to be roughly the scale at which gravity 
becomes perturbative, the universe could avoid a phase 
transition to the nonperturbative black hole phase, as¬ 
suming it begins in the weak-field phase. 


P(To,Ti) 



FIG. 2: Volume fraction in the weak-field phase of gravity 
as a function of the temperature To at which gravitons en¬ 
ter thermal equilibrinm. Shown are plots for three different 
valnes of the expansion speed, b = 10”^, 10”^, 10~® (left 
to right). Here the cut-off temperature Ti = 10“®, and the 
effective number of degrees of freedom g, = 10^. 


V. DISCUSSION 


We examined a possible first order phase transition 
of spacetime to a black hole phase with high entropy. 
Percolation of the high entropy phase occurs if gravi¬ 
tons are ever in a thermal state with temperature above 
Tg, either because they were born hot at the Planck 
epoch or because they were thermalized due to inter¬ 
actions with thermal matter. It seems possible, as sug¬ 
gested by entropic arguments, that almost all of gravi¬ 
tational phase space is accounted for by the nonpertur¬ 
bative phase. However, we find that the low entropy 
phase is metastable over timescales which are exponen¬ 
tially sensitive to the temperature, and potentially quite 
long. 

One may wonder how inflation changes this conclusion. 
We note that T§ is higher than the energy scale at which 
inflation is usually assumed to take place. If gravitons 
are never thermalized above a temperature of Tg then 
presumably inflation would simply take place as origi¬ 
nally envisioned. However, if gravitons are ever thermal 
with a temperature above Tg then we would speculate 
that it may be less probable for a given patch to inflate, 
although depending on the details of the model some non¬ 
zero probability may remain, even if a phase transition 
to the high entropy black hole phase does occur. 

Hot gravitons with temperature slightly below Tg will 
not lead to a phase transition; they will simply be red- 
shifted away. Both gravitons and matter may be born 
hot, as long as the temperature of the universe (either 
initially or after a period of inflation) is never greater 
than Tg. This does not require fine tuning because Tg is 
of the same order as (?pert, the energy scale below which 
quantum gravity effects are small. It may still be the 
case that the initial conditions represent a subset of mea¬ 
sure zero in the total phase space, which is dominated by 
the nonperturbative black hole phase |2fl|. However, our 
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analysis does show that once the initial choice of the low 
entropy phase is made, no transition to the high entropy 
phase need occur. 

These conclusions remain unchanged in a spacetime 
of arbitrary dimension d. For hot gravitons in d di¬ 
mensions, the exponent governing the nucleation rate 
of Eq. @ goes as F*/T ^ j'-{d- 2 ) ^ while for matter 
Flt/T ~ j'-(d-i){d- 2)/2 ^ Pqj. > 3 ^ black hole nucle¬ 
ation is suppressed more strongly in the matter system 
than in the gravitational one. Moreover, as d —> oo, Tg 
increases, meaning a transition to the black hole conden¬ 
sate phase is less likely in a universe with a large number 
of spacetime dimensions. 

Note added: After this paper was completed we be¬ 
came aware of earlier work using Euclidean path inte¬ 


gral methods in which Eqs. 0 , 0 , and 0 were inde¬ 
pendently derived j^. In these calculations imaginary 
time boundary conditions are applied to the gravitational 
field. Therefore, those authors were also studying ther¬ 
mal gravity and not only thermal matter. For discussion 
of black hole phase transitions, see 
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